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Abstract. We present a simple analytical description of the disruption of star clusters in a tidal field. The cluster 
disruption time, defined as tdis = {din M/dt} -1 , depends on the mass M of the cluster as tais = to(M / Mq) 1 
with 7 = 0.62 for clusters in a tidal field, as shown by empirical studies of cluster samples in different galaxies and 
by iV-body simulations. Using this simple description we derive an analytic expression for the way in which the 
mass of a cluster decreases with time due to stellar evolution and disruption. The result agrees excellently with 
those of detailed JV-body simulations for clusters in the tidal field of our galaxy. The analytic expression can be 
used to predict the mass and age histograms of surviving clusters for any cluster initial mass function and any 
cluster formation history. The method is applied to explain the age distribution of the open clusters in the solar 
neighbourhood within 600 pc, based on the new cluster sample of Kharchenko et al. that appears to be unbiased 
within a distance of about 1 kpc. From a comparison between the observed and predicted age distributions in the 
age range between 10 Myr to 3 Gyr we find the following results: (1) The disruption time of a 10 4 Mq cluster in 
the solar neighbourhood is about 1.3 ± 0.5 Gyr. This is a factor 5 shorter than derived from A-body simulations 
of clusters in the tidal field of the galaxy. Possible reasons for this discrepancy are discussed. (2) The present 
starformation rate in bound clusters within 600 pc from the Sun is 5.9 ± 0.8 10 2 MoMyr -1 , which corresponds 
to a surface star formation rate of bound clusters 5.2 ± 0.7 x 10 -10 Mq yr _1 pc~ 2 . (3) The age distribution of 
open clusters shows a bump between 0.26 and 0.6 Gyr when the cluster formation rate was 2.5 times higher than 
before and after. (4) The present star formation rate in bound clusters is about half as small as that derived from 
the study of embedded clusters. The difference suggests that about half of the clusters in the solar neighbourhood 
become unbound within about 10 Myr. (5) The most massive clusters within 600 pc had an initial mass of about 
3 10 4 Mq. This is in agreement with the statistically expected value based on a cluster initial mass function with 
a slope of -2, even if the physical upper mass limit for cluster formation is as high as 10 6 Mq. 

Key words. Galaxy: open clusters - Galaxy: solar neighbourhood - Galaxy: stellar content - Galaxies: star clusters 
- Stellar dynamics - 



1. Introduction 

Bound star clusters 1 in a tidal field are losing mass due 
to internal effects, i.e. mass loss by stellar evolution, and 
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Bound clusters are those that survive the infant mortality 
due to the removal of gas during the first 10 7 years (Fall 2004; 
Bastian et al. 2005) 
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by the external effect of tidal stripping. The combination 
of these effects results in a decreasing mass of the cluster 
until the cluster is destroyed completely. The time scale 
of this disruption depends on the initial conditions of the 
cluster, e.g. the stellar initial mass function and its concen- 
tration, and on the tidal forces experienced by the cluster 
during its galactic orbit. 

The disruption of star clusters determines the mass 
and age distributions of the existing clusters. Therefore 
any study of the cluster formation history of a galaxy has 
to take into account the disruption of clusters. 

The age and mass distribution of the initial cluster 
population is described by the cluster formation rate as 
a function of time, CFR(t), and the cluster initial mass 
function, CIMF. The distributions of the present observ- 
able star clusters is modified because: (a) the disruption 
time of the clusters depends on the initial mass, (b) the 
mass of each cluster decreases with time, (c) the clusters 
fade as they age due to stellar evolution. A (simple) de- 
scription of these three effects would facilitate the studies 
of samples of star clusters. The purpose of this paper is to 
provide such a simple description and show how it can be 
used in the analysis of star cluster samples. 

The structure of the paper is as follows: 
In Sect. 2 we discuss the arguments that the cluster dis- 
ruption time depends on its mass as a power law of the 
type idis oc M°- e2 . In Sect. 3 we describe the expressions 
used to calculate the evolution of a cluster in terms of its 
decreasing mass due to stellar evolution and disruption. In 
Sect. 4 we will show that the results agree very well with 
those of TV-body simulations of clusters in a tidal field. 
In Sect. 5 we use the description of the decreasing mass 
to predict the mass and age distributions of cluster sam- 
ples with various initial mass functions and various clus- 
ter formation histories. In Sect. 6 we apply the method by 
comparing the predicted age distribution of open clusters 
in the solar neighbourhood to the observed sample from 
Kharchenko et al. (2005). From this comparison we derive 
the disruption time of open clusters in the solar neigh- 
bourhood as well as the cluster formation rate and the 
star formation rate. We compare the results with inde- 
pendent determinations. The discussion and conclusions 
are in Sect. 7. 

2. A power law expression for the disruption time 
of star clusters 

Boutloukos & Lamers (2003, hereafter BL03) have stud- 
ied the mass and age distributions of magnitude limited 
cluster samples in selected regions in four galaxies, and 
concluded that these distributions can be explained if the 
disruption time of clusters depends on the initial mass Mi 
of the clusters as Mf , with 7 ~ 0.6 for clusters in very 
different local environments. 

Baumgardt & Makino (2003, hereafter BM03) have 
calculated a grid of N-hody simulations of clusters in the 
tidal field of our galaxy for different initial masses and 
initial concentration factors in circular and elliptical or- 



bits at various galactocentric distances. They take into 
account mass loss by stellar evolution and by tidal relax- 
ation. Their calculations show that the disruption time 
of a cluster, defined as the time when 5% of the initial 
number of stars remain in the cluster, scales with the half 
mass relaxation time i r h and the clusters crossing time 
i cr as tdis oc t* h tl~ x with x = 0.82 for clusters with an 
initial dimensionless depth Wo = 7 (which is a measure 
of the concentration index of the cluster, see King 1966) 
and x — 0.75 for less concentrated clusters with Wo = 5. 
BM03 and Giclcs et al. (2004) have shown that for all the 
models of BM03 the disruption time can be expressed as 
a function of the initial cluster mass as 

tdis = t (M,/M ) - 62 (1) 

where to is a constant that depends on the tidal field of the 
particular galaxy in which the cluster moves and on the 
ellipticity of its orbit. So the predicted dependence of the 
disruption time on the initial mass of a cluster agrees very 
well with the empirical relation derived by BL03. (De la 
Fuente Marcos & de la Fuente Marcos (2004) also report 
a power law dependence of the characteristic life time r 
of clusters on the number of stars. Based on a series of 
dynamical models they find that r <~ TV ' 68 where N is 
the initial number of stars of a cluster.) Lamers, Gieles & 
Portegies Zwart (2005) showed that to is expected to de- 
pend on the ambient density at the location of the clusters 
in that galaxy as to oc ' ■ 

The discussion above has concentrated on the compar- 
ison between the disruption time of clusters of different 
initial masses, i.e. t^ s oc M^ ' 62 , within one galactic envi- 
ronment. We have not yet discussed how the mass of an 
individual cluster decreases with time. This is the topic of 
the next section. 

3. The decrease of the cluster mass due to stellar 
evolution and tidal effects 

The mass of a cluster decreases due to stellar evolution 
and tidal disruption. We will describe the evolution of the 
bound mass, using analytic expressions for the mass loss 
from the cluster by stellar evolution and by tidal effects. 

3.1. Mass loss by stellar evolution 

The mass loss from clusters due to stellar evolution has 
been calculated for cluster evolution models by several 
groups, e.g. Bruzual & Chariot (1993) and the Starburst99 
models by Leitherer et al. (1999). We adopt the GALEV 
models for single stellar populations (Schulz et al. 2002; 
Anders & Fritze-v. Alvenslcben 2003). These models con- 
tain stars in the mass range of 0.15 < M* < 85 M@, 
distributed over this mass range with either the Salpeter 
or Scalo mass function. We adopt the models with the 
Salpeter mass function because deep photometry of clus- 
ters in the LMC shows that the cluster IMF is a power- 
law with a slope of about -2.35 down to at least 0.6 M Q 
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Table 1. Approximations to the mass lost by stellar evo- 
lution 

for GALEV cluster models with a Salpeter IMF of a = 
-2.35, 0.15 < M* < 85 M Q and 0.0004 < Z < 0.05 



z 


u-ev 






0.0004 


7.06 


0.265 


-1.790 


0.0040 


7.06 


0.260 


-1.800 


0.0080 


7.03 


0.260 


-1.800 


0.0200 


7.00 


0.255 


-1.805 


0.0500 


7.00 


0.250 


-1.820 



(de Marchi, 2003). Lower mass stars hardly contribute to 
the luminosity at ages less than 10 Gyr, but may con- 
tribute significantly to the cluster mass. The GALEV 
models are based on stellar evolution tracks from the 
Padova group, which include mass loss and overshooting 
(Bertelli et al. 1994; Girardi et al. 2000). Lamers (2005) 
has shown that the fraction of the initial cluster mass, Mj, 
that is lost by stellar evolution in the GALEV models, i.e 
q cv = (AM) CV /Mi where (AAf) cv is the mass lost by stel- 
lar evolution, can be approximated very accurately by a 
function of the form 

\ogq cv (t) = (log* - a ov ) b " + c cv for t > 12.5 Myr (2) 

The values of a cv , b cv and c cv are listed in Tabled for dif- 
ferent metallicities. This function describes the mass loss 
fraction of the models at t > 12.5 Myr with an accuracy of 
a few percent. The mass loss at younger ages is negligible 
because the most massive stars with M» > 3OM hardly 
contribute to the mass of the cluster. For cluster models 
with a lower limit of the stellar IMF different from 0.15 
M Q , the value of q cv (t) can easily be adjusted, because 
stars with M < 0.6M Q contribute to the cluster mass 
but not to its mass loss at ages less than 10 10 yrs. (This 
mass loss rate is very different from that of the Starburst99 
models, because the Starburst99 models have a lower limit 
for stellar mass of 1 Mq.) 

In this paper we use the symbol /j,(t) = M{t)/Mi to 
describe the fraction of the mass of a cluster with initial 
mass Mi that is still bound at age t. We define 



/i ev (t) = 1 - q ev (t) 



(3) 



as the fraction of the initial mass of the cluster that would 
have remained at age t, if stellar evolution would have 
been the only mass loss mechanism. The function (i cv (t) 
is independent of the initial mass of the cluster. 

3.2. Mass loss by stellar evolution and tidal effects 

We describe the decreasing mass of a bound cluster that 
survived infant mortality (t <; 10 7 yrs) as a function of 
time. Let us define M(t; M\) as the mass of a cluster of 
initial mass Mi and age t, and n(t; Mj) = M(t; M\)jMi as 




Fig. 1. The predicted decrease in cluster mass due to stel- 
lar evolution and disruption for Z — 0.020 and four values 
of the initial cluster masses: 10 3 , 10 4 , 10 5 and 10 6 Mq. 
Top panel: to — 2 Myr; lower panel: to — 30 Myr. The 
full lines give the exact decrease derived by numerical so- 
lution of the differential equation (@J , and the dotted lines 
give the approximation (Eq.EJ. Notice the excellent agree- 
ment. The dashed lines gives the decrease in mass due to 
stellar evolution only. 

the fraction of the initial mass that is still in the cluster. 
The decrease of mass due to both stellar evolution and 
disruption can then be described as 



dM _ fdM\ (dM^ 



(4) 



where the first term describes the evolution by stellar mass 
loss and the second term by disruption. Following the ar- 
guments given in Sect. 2 we assume that we can describe 
the mass loss by disruption as 



(dM 
\dt 



-M 



-M 



-Mr, 



dis 



t dis f (M/M )7 t c 



with 7 = 0.62 and to is a constant that depends on 
the tidal field. (See Lamers, Gieles & Portegies Zwart, 
2005 for the dependence of to on the conditions in dif- 
ferent galaxies). The first equality assumes that the mass 
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lost by disruption can be approximated by an exponen- 
tial decay with a time scale that decreases as the mass of 
the cluster decreases. This is equivalent to the statement 
that the disruption time in our description is defined as 



tdit 



dlriM/dt. The second equality assumes that this 



timescale depends on the mass as M 7 . (If 7 was equal to 
1 and the evolutionary mass loss could be ignored, then 
the mass of the cluster would decrease linearly with time 
until t = toMi/M®. For 7 = the decrease would be 
exponential.) 

Equation 0] can easily be solved numerically. It turns 
out that the mass decrease of a cluster can be approxi- 
mated very accurately by the following formula 



M(t) 



< 11 ( M ® 



71 1 h 



(6) 



if the first term in brackets is larger than the second 
term. If the second term is larger than the first term, 
i.e. when the mass lost by disruption is larger than the 
mass that remained after mass loss by stellar evolution, 
then fi(t; Mi) = and the cluster is completely disrupted. 
Approximation^^ quite accurate because during the first 
10 8 years mass loss is dominated by stellar evolution so the 
second term is negligible and fJ, cv (t) describes the fraction 
of the mass that survives mass loss by stellar evolution. 
During later years, when /i ev (t) decreases very slowly, the 
mass loss is dominated by disruption, which is described 
by the second term in brackets. Eq. Elcan be inverted to 
express the initial cluster mass in terms of the present 
cluster mass: 



Mi 



M 
~Mq 



to 



1/7 



(7) 



Figure ^ compares the numerical solution with the 
analytic approximation for various initial masses, 10 3 < 
Mj/M© < 10 6 , for a short and a long disruption timescale, 
t = 2 Myr and 30 Myr, both for Z=0.02. In all cases the 
agreement between the analytic and the numerical solu- 
tion is excellent, i.e. within about 0.015 dex, although the 
disruption times vary by more than 6 orders of magnitudes 
from model to model. Even for the low mass cluster model 
of Mi = 10 3 M© at to — 2 Myr, for which the disruption 
is already effective during the first 10 Myr, the agreement 
between the the numerical solution and the analytic ex- 
pression is very good. (Tests show that Eq.[S]is also a very 
good approximation for all other values of 7 in the range 
of < 7 < 1.) 

We define t^ al as the total disruption time and 
as the time when only a fraction A of the initial mass 
remains. ^From Eq. we find that ^° s tal and ^r A are 
described by the implicit relations 



(lis 



and 



l dis 



total _ to f Mj 1 , . totals 17 

~ ~ 1 M© I i/iev ^ dis n 



Mi 

M Q 



A 7 ) 



(8) 



(9) 



We will use this last expression for a comparison of our 
analytic solution with those of ./V-body-simulations. In the 
range of 10 4 < t < 10 7 years and 10 3 < M z < 10 6 M© 
can be approximated by 



the values of t^°* al 



log^total) 



log f^-J +7I0 
-0.00825 log 



Mj 
M Q 



Mj 
Mq 



x log 



10 4 yr 



(10) 



This approximation is valid within 0.03 dex for all metal- 



licities. For all models we find that t^° 



0.89 i£* al . 



These equations imply that the total disruption time 
of a cluster with an initial mass of 10 4 M©, which was 
referred to as t\ in BL03 and in Lamers, Gieles & Portegies 
Zwart (2005), is related to to by 



total 



1.355 IO 47 



-t. 



0.967 



7 



6.60 x 10 2 t 



0.967 



(11) 



where the last equation is only valid if 7 = 0.62 and to is 
in yrs. 

Equation El shows that t^°* al (Mi) is approximately 
proportional to M^ 7 , which was the relation adopted in 
the study of the disruption times of clusters in different 
galaxies by BL03 and Lamers, Gieles & Portegies Zwart 
(2005). 

4. Comparison of the analytic solution with 
results of iV-body simulations 

We can compare our analytic expression for the decreasing 
mass of a cluster with the results of TV-body simulations. 
We adopt the simulations by BM03 who calculated the 
fate of clusters under various conditions. Before making 
this comparison, we would like to point out that: 

(a) BM03 adopted the stellar IMF of Kroupa (2001), 
which have an initial mean mass of 0.547 Mq, whereas our 
calculations are based on the GALEV cluster evolution 
models (see Sect. 3.1) which have a Salpeter (1955) IMF 
with a lower mass cut-off of 0.15 M© and a maximum 
mass of 85 M©, resulting in an initial mean mass of 0.516 

M©. 

(b) BM03 adopted the evolutionary mass loss rates 
from Hurley, Pols & Tout (2002), whereas those of the 
GALEV models are based on the calculations from the 
Padova-group (see Sect. 3.1). 

(c) BM03 define the disruption time of a clusters, 
(called dissolution time in their paper) as the time at 
which only 5 percent of the initial mass is still in the 
cluster. 

(d) The simulations by BM03 show that the mean mass 
of the remaining stars in a dissolving cluster changes with 
time. In the early phase the mean mass decreases because 
stellar evolution removes the massive stars, but the mean 
mass increases in later phases when disruption becomes 
the dominant mass loss mechanism and low mass stars 
are lost preferentially. In the simulations by BM03 the 
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-0.6 -0.4 -0.2 0.0 0.2 -0.6 -0.4 -0.2 -0.0 0.2 

l°g 'Am* log t/t^s 

Fig. 2. Comparison of the decrease of the cluster mass with time between the results of the iV-body simulations by 
BM03 (left) and our description (right) for clusters of different initial numbers of stars, Ni or mass M- x . The mass 
has been corrected to first order for the mass lost by stellar evolution. In all figures the time t is scaled to t^°, i.e. 
the time when 95% of the cluster mass is lost due to stellar evolution and disruption. The upper figures are for mass 
versus time and the lower figures are for log(Mass) versus log(time). Notice the strong similarity between the results 
of the ./V-body simulations and our simple description (Eq. 6). 



mean mass near the end of the life of a cluster of initially 
3 x 10 5 stars has increased from 0.516 to 1.2 M Q . In 
our analytic approximation this effect is not taken into 
account Therefore we can expect a small offset in the 
timescale between the iV-body and the analytic results 
due to mass segregation. 

Because of differences in the adopted mass loss by stel- 
lar evolution between our and BM03 models, we compare 
the results for the decreasing mass of a cluster not di- 
rectly, but corrected for the stellar evolution. This means 
that we compare the predictions for ^i(t; Mi)/ /x ev (i) rather 
than for n(t;Mi). The function /j,(t; Mi) / /j, ev (t) describes 
the fraction of the initial mass of the cluster that is lost 
by disruption only. This function, which is called M re \(t) 
by BM03, is expected to be approximately independent of 
the adopted evolutionary mass loss rates. 



BM03 give the function M re \(t) for their models of clus- 
ters of different initial numbers of stars, 8.2 x 10 3 < N[ < 
1.3 x 10 5 , and with an initial concentration described by 
Wq = 7, in circular orbits around the galactic center at a 
distance of 8.5 kpc. (sec BM03 Fig. 6). Their results show 
that M rc \(t) decreases almost linearly with time, and that 
the function plotted against t/t^^f is about the same for 
all clusters. We compare their results with our analytic 
expression Eq. El 

The left hand panels of Fig. [2] show the normalized 
decrease in cluster mass, M IC \(t / t^^f) , of the models by 
BM03. The top figure gives the mass as function of time, 
both in linear scale. However, since the mass of clusters 
will decrease orders of magnitudes before they are dis- 
rupted, we also plot the logarithm of M IC \(t) as a func- 
tion of the logarithm of t/t^f in the lower figure. The 
right hand panels of Fig. [21 show the results of our mod- 
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els. For a fair comparison between our calculations and 
those of BM03, we plot fj,(t;Mi)/fj, ev {t) as a function of 
^Adfs % - The agreement between the predictions of the N- 
body calculations and our description is very good, both 
in linear scale as in logarithmic scale. The very small dif- 
ference of less than 5 percent during the early phase of 
the most massive models is probably due to the difference 
in the adopted mass loss by stellar evolution. (The plot- 
ted relations are to first order corrected for the effects of 
stellar evolution. However because the cluster evolution is 
caused by the combination of stellar evolution and tidal 
effects, the corrected relations still bear the imprint of 
the adopted stellar evolution.) Figure |3 shows the direct 
comparison between the A^-body prediction for a clusters 
of initial mass M, = 8.9 x 1O 3 M (N t = 16384) with a 
concentration of Wq — 7 and our analytic solution. For 
the analytic solution we adopted the same mass and the 
parameter to = 18 Myr was chosen in such a way that 
the 95 % disruption time is 5.7 Gyr, very close to that of 
the iV-body simulation. The times are normalized to the 
time when 95% of the clustermass is gone. We see that the 
prediction by the TV-body calculations (full line) and the 
analytic expression (dotted line) are very similar, apart 
from a small offset of the timescale. If we normalize the 
timescale of the analytic solution to the time when 96.5% 
of the initial mass is gone (dashed line), the agreement is 
almost perfect! This difference in timescale is due to the 
fact that in the TV-body simulations there is a preference 
for the low mass stars to be kicked out of the clusters, 
whereas in the analytical solution stars of all masses are 
lost. Apart from this small difference, the analytical solu- 
tion with to as a free parameter describes the decrease of 
the mass of clusters surprisingly accurately. 
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Fig. 3. Comparison between the mass decrease of a cluster 
of 8.9 x 1O 3 M predicted by N-body simulation of BM03 
(full line) and our analytic approximation for to — 18 Myr 
(dashed and dotted line). The mass decrease has been cor- 
rected for the mass lost by stellar evolution. The timescale 
is normalized to the time when 95 or 96.5% of the initial 
mass is lost (see text). 

ing time t CI as about tf h tl~ x , with x ~ 0.75 for models 
with a concentration parameter Wq = 5.0 and x — 0.82 
if W = 7.0 (BM03). For both sets of models the disrup- 
tion time scales with mass as M 062 (Gieles et al. 2004). It 
is the continuous adjustment of the half mass relaxation 
time and the crossing time to the changing conditions of 
the cluster that results in an exponential disruption time 
that varies during the life of a cluster as the present mass 

M 0.62 



4.1. Explanation 

We conclude that our analytic description of the decreas- 
ing mass of a cluster, with the adjustable free parameter 
to, is very similar to the one derived by TV-body simula- 
tions. We want to stress that this result is not trivial. 

The values of t^f that results from the A^-body simu- 
lations by BM03 depend on the initial number of stars as 
A~j°- 62 , in excellent agreement with the empirically derived 
dependence by BL03 (see Sect. 2). However this does not 
automatically imply that the decrease of mass with time 
of each individual cluster can be described by a function 
that depends on the present mass of that cluster to the 
same power 7 = 0.62 (Eq.EJ. However the good agreement 
between the simulations by BM03 and our result from the 
analytic description shows that it does. 

So, not only does the total lifetime of all clusters de- 
pend on their initial mass as t*^ 1 oc M^ ' 62 , but also at 
every moment during the lifetime of a cluster the exponen- 
tial disruption time, (din M/dt)" 1 , scales with the mass 
to the same power 0.62. This is a consequence of the fact 
that the timescale of the disruption process depends on 
both the half mass relaxation time, i r h and the cross- 



5. The predicted mass and age distributions of 
cluster samples 

Using the expression for the decreasing mass of clusters, 
Eq. El we can predict the mass and age distribution of 
cluster samples (for open clusters as well as globular clus- 
ters) for any adopted cluster formation rate, CFR(t), and 
cluster initial mass function, CIMF. 

Suppose that the CIMF is a power law with a slope 
-a = -2 (Zhang & Fall 1999; Larsen 2002; Bik et al. 2003; 
de Grijs et al. 2003) in the range of M m - m < M c \ < M max , 
TWmin ~ 1O 2 M and M max rj 10 7 M q , then the number of 
clusters with initial mass Mj formed at time t will be 

N(M u t) = S(t) (j^-j for M min < Mi < M max (12) 

in Nr M _1 yr~ 1 if t is in years and M, in M . The func- 
tion S(t) is related to the cluster formation rate CFR(t) 
in Nr yr -1 as 

/■TWmax 

CFR(t) = / N(M h t)dMi 
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S(t) 



'Mr, 



Mr. 



M„ 



Mr. 



(13) 

The total mass of the clusters formed per year is 
S(t) ln(M max /M min ) in M Q yr" 1 for a = 2. 

5.1. The distribution of the masses and ages 

Given the initial mass distribution of the clusters, their 
formation rate, CFR(t), and an expression for the way 
in which the mass of each cluster changes with time (Eq. 
EJ), we can calculate the present distribution of existing 
clusters as a function of age or as a function of their mass. 

If N(M,t) is the number of clusters of mass M and 
age t in Nr M© -1 yr _1 , then N(M,t) and N{M u t) are 
related by the conservation of the numbers of clusters 



N(M, t) dM = N(Mi,t) dMi 



(14) 



with M(t) and Mi(t) related via Eq. El Applying the 
derivative d(M,t)/d(Mi,t), that follows from Eq. and 
combining this with Ea.[Plfor the CIMF and Ea.lT31for 
the CFR we find the present distribution of clusters as 
function of mass and age: 



N(M,t) = S(t) 



' M 






r 


It ( 


M 


to V 


~Mq 



A*ev(i) 



a-l 



-i)h 



(15) 



This equation is valid for M smaller than some upper 
limit, M up (t), which is the mass of a cluster of age t with 
the maximum initial mass M lrl! , x 



M up (t) = M n 



to 



Mq 
M ma; 



1/7 



(16) 



Similarly, for a given value of M, Eq. 1151 is only valid for 
ages less than t up (M) which is the age at which a cluster 
with an initial mass of M max has reached a mass M. So 
t up (M) is given by the condition 



M 



to 



M r . 



(p cv (tu P )V = (17) 



Eq. E| allows us to calculate the predicted mass and age 
distribution of a cluster sample for any assumed cluster 
formation rate. The mass distribution is found by inte- 
grating N(M, t) over age for any mass, and the age distri- 
bution is found by integrating over mass between M up (t) 
and some lower mass limit Mi ow (t), set by the detection 
limit, for any age. 

The mass and age distribution of cluster samples 
N(M, t) depends on the stellar evolution and on disrup- 
tion. To get insight into the effect of disruption on the 



21 
~25~ 



log N(M,t) 

no Mdot(evol): // ev =l. 
parameters: log t/t 

no disr : slope = - a 
strong disr: slope = y— 1 




log (M/M e ) 

Fig. 4. The changes in the mass distribution (Eq. 15) of a 
sample of clusters as a function of their age, in case stellar 
evolution can be neglected. We adopted a cluster initial 
mass function in the range of 10 2 < M < 10 6 Mq with 
a — 2.0 and a disruption parameter 7 = 0.62. The differ- 
ent curves refer to different ages, which are parametrized 
by log t/to- The maximum mass decreases with age due 
to disruption. 



log Ni„,/S 

no Mdot(evol): A i e v = l- 
M m „=10 6 M Q 

parameters: log M mjn /M 




log t/t 

Fig. 5. The age distribution of cluster samples formed at 
a constant formation rate, in case stellar evolution can be 
neglected. We adopted a cluster initial mass function in 
the range of M min < M < M max with M max = 10 6 M Q 
and different values of M m in- We adopted a = 2.0 and 
a disruption parameter 7 = 0.62. The curves are labeled 
with log (M min /M Q ). 

evolution of a cluster sample we first consider a simpli- 
fied case when mass loss by stellar evolution is neglected, 
i.e. fi cv (t) = 1.0. In that case the function N(M,t)/S(t) 
depends only on the slope a of the cluster IMF, the mass- 
dependence 7 of the disruption and on the ratio t/to- 
FigureHshows the shape of N((M, t) /S(t). For very young 
ages or very long disruption time (t/to < 10) the distri- 
bution is the initial CIMF with slope —a. For strong dis- 
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ruption, i.e. t/to > 10 2 , the distribution of the low mass 
clusters becomes flatter and approaches a power law of the 
type N(M) ~ A/ 7-1 . This distribution is similar to the 
one predicted by BL03 for instantaneous disruption, ex- 
cept that the transition between the two slopes is gradual, 
whereas it shows a sharp kink for models with instanta- ; 
neous disruption. 

Figure |S] shows the normalized age distribution 
Ntot(t)/S, for cluster samples formed at constant cluster 
formation rate in the mass range of M m j n < M < Af max 
for various values of M m i n , in case the mass loss by stellar 
evolution can be ignored. The distribution is flat for young 
clusters at a value of N/S ~ M~ ia an d curves down to 
older clusters, approaching a slope N ~ (t/t )~ 1 /' 1 . This 
was predicted by BL03 for instantaneous disruption. The 
distribution drops to zero at the age at which the most 
massive clusters are disrupted, i.e. when t/to = ^maxll 
(Eq. 16), which is at t/to = 8.46 10 3 for M max = 1O 6 M 
and 7 = 0.62. 

6. Application to Galactic open clusters 

Ideally one would like to compare the predictions with 
complete (or at least unbiased) samples of clusters with 
known masses and ages. Unfortunately this is not possi- 
ble at the moment, because samples of clusters in external 
galaxies are usually magnitude limited. (The method for 
determining the disruption times from magnitude limited 
cluster samples with gradual disruption will be described 
by Lamers (2005) and applied to the cluster sample in 
M51 by Gieles et al. (2005a)). Samples of open clusters in 
the solar neighbourhood are unbiased, but only the clus- 
ter ages have been determined systematically and not the 
cluster masses. We will compare our predictions for the 
age distribution to the sample of open clusters in the so- 
lar neighbourhood. 

6.1. The sample of open clusters 

Kharchenko et al. (2005) have recently published a cat- 
alogue of astrophysical data of 520 galactic open clus- 
ters (COCD = Catalogue of Open Cluster Data) in the 
wider neighbourhood of the Sun with the values of angular 
sizes of cluster cores and coronae, heliocentric distances d, 
E(B—V), mean proper motions, radial velocities and ages. 
These parameters have been determined by homogeneous 
methods and algorithms including a careful procedure of 
cluster member selection. The basis of this study is the 
ASCC-2.5 - All-Sky Compiled Catalogue of about 2.5 mil- 
lion stars (Kharchenko 2001) down to V ~ 14 (complete- 
ness limit at V ~ 11.5), with compiled proper motions and 
B, V magnitudes, which is based on the Tycho — 2 data, 
and supplemented with Hipparcos data sets, as well as 
with some ground-based catalogues 2 . Cluster member- 
ship is based on a combined probability which takes into 




2.5 3.0 
log(d) [pc] 

Fig. 6. Top: the surface density distribution projected 
onto the Galactic plane of open clusters in the solar neigh- 
bourhood from the homogeneous Kharchenko et al. (2005) 
catalogue. Errorbars indicate la statistical uncertainties. 
The surface density is about constant up to at least 600 pc 
and possibly 1 kpc. Bottom: the ratio between the num- 
bers of old (> 2.5 10 8 yr) and young (< 2.5 10 8 yr) clusters 
as a function of distance. The ratio is about constant up 
to 1 kpc. 



Milky Way d<600 pc 
1 14 clusters 

Kharchenko ct al. (2005) 



2 The ASCC-2.5 catalogue can be retrieved from the CDS at 
|ftp:/ / cdsarc.u-stra sbg.fr / pub /cats/I/280 A 



log Age [yr] 

Fig. 7. The age histogram in units of number per year, 
in logarithmic age-bins of 0.2 dex, of 119 open clusters 
within d < 600 pc from Kharchenko et al. (2005). The 
distributions are plotted for two sets of bins, shifted by 
0.1 dex, with and without squares respectively. The error- 
bars indicate the la statistical uncertainty. The distribu- 
tion decreases to older ages, with a small bump around 
log(t/yr) ~ 8.6. For log(f/yr) < 7.5 the distribution is 
uncertain due to large errorbars. 
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account kinematic (proper motion), photometric, and spa- 
tial selection criteria (see Kharchenko et al. 2004 for de- 
tails). For stars within a circle with a cluster radius the 
membership probability is calculated as a measure of a de- 
viation either from the cluster mean proper motion (kine- 
matical probability), or from the Main Sequence edges 
(photometric probability). Stars deviating from the refer- 
ence values by less than one a (rms) are classified as most 
probable cluster members (la-members, i.e., with a mem- 
bership probability P > 61%). Those falling in semi inter- 
vals [lcr,2(7) or [2cr, 3cr) are considered as possible members 
(P = 14 - 61%) or possible field stars (P = 1 - 14%), 
respectively. Stars with deviations larger than 3cr are re- 
garded as definite field stars (P < 1%). As a rule, all clus- 
ter parameters were determined from the data of the most 
probable cluster members. Cluster ages were determined 
with an isochrone-based procedure which provides a uni- 
form age scale (see Kharchenko et al. 2005 for details). 
Thus the COCD is the most homogeneous and most com- 
plete catalogue of open clusters in the solar neighbourhood 
available up to now. 

Figure shows the distance distribution of the density 
of clusters projected onto the Galactic plane, in number 
per pc 2 . We see that within the statistical uncertainty the 
surface density is constant up to at least 600 pc, and possi- 
bly even up to 1 kpc. The lower part of the figure shows the 
ratio between old (t > 2.5 10 8 yr) and young (t < 2.5 10 s 
yr) clusters as a function of distance. Up to a distance 
of about 1 kpc there is no significant change in this ratio 
within the statistical uncertainty. This is important for our 
study, because it shows that the age distribution of open 
clusters within about 1 kpc is not affected by detection 
limits. 

Figure[7]shows the age distribution in number per year 
of the 114 clusters within 600 pc in the Kharchenko et 
al. (2005) sample. The effect of binning is demonstrated 
by plotting two sets of data, where the bins have been 
shifted by 0.1 dex relative to one another. This distribu- 
tion is decreasing with age, apart from a small local max- 
imum around log(f/yr) ~ 8.5. The distribution at young 
ages is sensitive to the choice of the age-bins and shows 
a significant scatter. The steep slope at log(i/yr) > 8.8 
demonstrates that cluster disruption is important. 

6.2. The lower mass limit of the clusters 

For the determination of the cluster formation rate and 
the disruption times we need an estimate of the minimum 
mass of the clusters in the Kharchenko et al. (2005) sam- 
ple. This catalogue does not list the mass of the clusters, 
but it can be estimated roughly from the age, distance, 
extinction and the number of stars of each cluster. We 
have estimated the lower mass limit of the Kharchenko et 
al. cluster sample in the following way. 

(a) First we calculate the number of members brighter 
than the completeness limit, Vy lm — 11.5, within the clus- 
ter radius with the following constraints on cluster mem- 



bership probabilities: 2a photometric probability and 2a 
kinematic probability. These probabilities were defined in 
Sect 6.1. We did this separately for main sequence stars 
only, and for members of all spectral types. 

(b) Using the distance and E(B — V) of the clusters, we 
expressed Vj; m in M v . This limiting magnitude of the main 
sequence (MS) stars is expressed in M^ \ and mass, M* im , 
using the bolometric corrections and the mass luminosity 
relation of luminosity class V stars. With the cluster age 
known, the mass of the stars at the turn-off point of the 
MS, MtO) can be estimated from the relation between the 
MS lifetime and stellar mass, for which we adopted the re- 
lation by Schaller et al. (1992) for solar metallicity. 

(c) We then assumed a stellar IMF with a slope of -2.35, 
i.e. N(M)dM = CM~ 2 - 35 dM, and calculated the value of 
C that gives the derived number of main sequence stars 
in the mass range of M* im < M < Mro- 

(d) With this value of C we calculated the total mass of 
the cluster for all stars between the upper MS mass limit, 
Mto, and a lower limit for the stellar mass, M^ in , for 
which we adopted 0.15M©. (With this lower mass limit 
the mean stellar mass of a cluster with a Salpeter mass 
function is 0.51Af Q , which is quite similar to the mean 
mass of 0.55M© for a Kroupa (2001) mass function.) We 
corrected this mass for the small number of stars with a 
MS age that is 15% shorter than the age of the cluster, in 
order to correct for the stars that have evolved off the MS 
but have not yet ended their live. So the adopted mass 
range is M^c < M < ^mim where M^, is the mass of 
a star with a MS lifetime of 0.85 times the age of the clus- 
ter. (We note that white dwarfs, neutron stars and black 
holes do not add significantly to the mass of clusters with 
ages less than about a few Gyr.) In this way we estimated 
the mass of all the clusters in the Kharchenko sample with 
d < 600 pc. 

(e) We also applied this method directly to the observed 
number of the probable (2a) member stars with V < 
11.50, of all luminosity class. These are the observed stars 
in the mass range of < M < M* lin . The resulting 
masses are very similar to those derived from the number 
of MS stars only, except for a few clusters of high extinc- 
tion for which the Vu m — 11.5 corresponds to stars near 
the top of the main sequence and the number of probable 
members brighter than V — 11.5 is small. 

(f) To estimate the sensitivity of the resulting cluster mass 
to the adopted stellar lower mass limit, we repeated the 
analysis for an adopted lower mass of M* nin — 0.25M©. 
In this case the estimated masses are about 80% of those 
estimated for M* ln = 0.15M©. 

The resulting mass-age histogram of the clusters, for 
M^ in = 0.15M©, is shown in Fig. El Most of the clusters 
have a present mass in the range of about 5 10 1 to 5 10 3 
Mq. The histogram of the resulting masses, Fig. shows 
a peak in the range of about 100 to 300M©. The slow 
decline to the high mass end reflects the initial cluster 
mass function modified by mass loss and disruption. The 
steep decrease to low masses is due to the detection limit 
of the clusters and their members. The edge suggests that 



10 



H.J.G.L.M. Lamers et al.: Disruption of star clusters in tidal fields 



□ Dm o D S 

□ □ r. 



□ 



log Age [yr] 

Fig. 8. The mass-versus-age diagram of 114 clusters of 
the Kharchenko et al. (2005) catalogue within a distance 
of 600 pc. The mass is derived from the number of main 
sequence stars with V < 11.50. 



Kharchenko sample 



d < 600 pc 
114 clusters 



log Mass [M G ] 

Fig. 9. The mass histogram of 114 clusters of the 
Kharchenko et al. (2005) catalogue within a distance of 
600 pc. The steep edge at the low mass side suggests 
that the sample is complete for clusters with a mass 
M k, WOMq. 



the mean lower mass limit of the Kharchenko et al. (2005) 
cluster sample is about IOOA/q. This mass corresponds 
to a minimum number of about 280 stars per cluster if 
M* nin — 0.15M Q . From Eq. 0we find that a present mass 
of 100 Mq corresponds to an initial mass of 3.4 10 2 Mq if 
t = 10 8 yr and 5.8 10 3 M if t = 10 9 yr. These values are 
for a disruption parameter of t$= 3.3 Myr (see below). 

To estimate the sensitivity of the cluster masses to 
the adopted stellar lower mass limit we also determined 
the masses of the clusters in the Kharchenko et al. (2005) 
catalogue within 600 pc in the same way as described 
above but with an adopted minimum stellar mass of 
M^ in = O.25M instead of 0.15 M Q . The resulting cluster 



masses are about 80% of those for M^ in = 0.15M Q , so in 
that case the minimum lower mass limit of the Kharchenko 
sample would be about 80 Mq. The limiting cluster mass 
of 80 Mq corresponds to about 140 stars. 



6.3. The disruption time of clusters in the solar 
neighbourhood 

Figure HHI shows the fits to the data for 7 = 0.62 and var- 
ious values of to, based on the method described in Sect. 
15.11 For the top figure we adopted a constant cluster for- 
mation rate and a CIMF with a — 2 and a mass upper 
limit of 1 10 5 Mq. This latter choice agrees with the steep 
decrease around log(t/yr) ~ 9.5. We assumed a minimum 
detectable cluster mass of 100A/ Q . The predicted distribu- 
tions are normalized to the data point at log(i/yr) = 8.1, 
which is one of the most accurate data points. The best fit 
is reached for to — 3.3 Myr. (We remind that the low dat- 
apoint at log(i/yr) = 7.3 and the subsequent high point 
at 7.5 may be due to the adopted binning; see Fig. 0) 

The distributions for to >> 3.3 Myr underpredict the 
observed numbers at log(f/yr) < 7.2 and overpredicts the 
numbers at old ages. The distribution for to << 1.6 Myr 
overpredicts the distribution at young ages. We have ap- 
plied a x 2 test to express the goodness of the fit. For this 
test we only considered the data at 7.1 < log(i/yr) < 8.1 
and log(i/yr) > 8.9, i.e. we excluded the bump at 8.3 < 
log(t/yr) < 8.7 which will be discussed below. We also ex- 
cluded the data younger than 10 Myr, because they may 
be affected by infant mortality. The criterion \ 2 < Xmin + 1 
results in a best estimate of to = 3.1^J' § Myr. The data 
of the fits are given in the top line of Table [5] 

To investigate the effect of the adopted lower mass 
limit, we have also compared the observed age distribution 
with that predicted for M max = 1 10 5 Mq and M min = 80 
Mq . This last value is derived from the Kharchenko et al. 
cluster sample if the stellar lower mass limit of 0.25Mq 
is adopted instead of 0.15Mq. This comparison is very 
similar to the one in the top panel of Fig. and is not 
shown here. The resulting data are listed in the second line 
of Table [5] The best fit is at a slightly longer disruption 
time of to = 3.5*q g Myr. Combining the two values of 
to, derived for M min =80 and 100 M , we conclude that 
t = 3.3^;;; Myr (see bottom line of Table©. 

The derived value of to implies a total disruption time 
of a 10 4 Mq cluster of 1.3 ±0.5 Gyr (Eq.ETJ. This empir- 
ically derived disruption time of open clusters is about a 
factor 5 smaller than the value predicted by the iV-body 
simulations of BM03, who predict a disruption time of 6.3 
Gyr for a 10 4 Mq cluster with an initial concentration 
factor of Wo = 5.0 and 5.9 Gyr if W = 7.0. The iV-body 
simulations of clusters by Portegies Zwart, Hut & Makino 
(1998) resulted in about twice as short disruption times as 
those of BM03 (see also Lamers, Gieles & Portegies Zwart, 
2005) . This is still longer than the empirically derived dis- 
ruption time for open clusters in the solar neighbourhood. 
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Fig. 10. Comparison between observed and predicted age 
histogram of clusters in the solar neighbourhood within 
d < 600 pc. The data are fitted to predicted relations 
based on our analytical expression of the cluster disrup- 
tion with various values of to , normalized to the point at 
log(i/yr) = 8.1. The clusters are formed in the mass range 
of 10 2 < M cl < 10 5 Mq, with a CIMF of slope -2.0. Top 
figure: Predictions for a constant CFR. The dotted line 
indicates the prediction if there was no cluster disruption, 
but only mass loss by stellar evolution. The shorter the dis- 
ruption time, the steeper the decrease towards high ages. 
At young ages the shortest disruption time corresponds to 
the largest formation rate and vice versa. Lower figure: the 
best fit for an assumed burst between 250 and 600 Myr 
ago, with a CFR that was 2.5 times higher than before 
and after the burst. 

We consider three possible reasons for this discrepancy: 
(a) the clusters do not start with the initial concentration 
factors adopted in the simulations by BM03 and (b) the 
presence of another mechanism (apart from the tidal field) 
that contributes to the destruction of clusters in the solar 
neighbourhood. 

BM03 assumed in their iV-body simulations that the 
clusters initially fill their tidal radius with a concentration 



Table 2. The disruption time and the formation rate of 
clusters within 600 pc from the Sun for two assumed values 
of the lower mass limit of clusters in the Kharchenko et 
al. (2005) sample. 



Mass range 
(M ) 


to 
(Myr) 


log(CFR) 
Nr/yr 


log(CF7?) 
M /yr 


8 10 1 - 1 10 5 


o tC+1-3 

°- d -i.o 


-6.03 ± 0.02 


-3.27 ±0.02 


1 10 2 - 1 10 5 


o -,+1.2 
J -- L -0.8 


-6.03 ± 0.02 


-3.19 ±0.02 


adopted 


o q+1.4 
o - J -1.0 


-6.03 ± 0.02 


-3.23 ± 0.06 



factor of Wo = 5.0 or 7.0, defined by King (1966). These 
are the values suggested by the current density profiles 
of the globular clusters. However, open clusters are much 
less centrally condensed than globular clusters. Moreover 
they may not fill their tidal radius when they are formed. 
If the clusters are smaller than their tidal radius the in- 
ternal relaxation will be faster than predicted and so the 
disruption might be faster than predicted. We suggest that 
open clusters are formed so far out of equilibrium that they 
lose a substantial fraction of their mass within a few cross- 
ing times. Most clusters will then disperse completely, in 
agreement with the high infant mortality rate. The surviv- 
ing clusters might then dissolve along the lines predicted 
by the A-body simulations, but on a faster time scale. 
A-body simulations of clusters with various initial con- 
centration factors and various initial radii arc needed to 
test this suggestion. 

The disruption times calculated by BM03 is an 
upper limit because the the values are calculated for 
tidal disruption in a smooth tidal field without other 
destruction mechanisms. The destruction of open clusters 
by encounters with giant molecular clouds (GMCs) has 
been proposed by several authors, e.g. Oort (1958) and 
Terlevich (1987). The problem with this explanation is 
that disruption by a CMC is expected to result in a 
mass dependence disruption of tdis ~ M 1 with 7 = 1.0 
(e.g. Spitzer 1987), whereas the observed age distribution 
agrees perfectly with 7 = 0.62 predicted for tidal effects. 
On the other hand, massive clusters are likely to have 
larger tidal radii and hence more interactions than low 
mass clusters and will therefore be more susceptible 
to the influence of GMCs. This might soften the mass 
dependence of the disruption time to 7 < 1. 



6.4. The cluster formation rate, the star formation rate 
and the infant mortality rate 

The vertical shift of the predicted relative to the observed 
age distribution in the top panel of Fig. HUI gives the CFR. 
We find a CFR of 0.93 ± 0.04 clusters per Myr. (This is 
about twice as high as the value derived by Battinelli & 
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Capuzzo-Dolcetta (1991) based on the Lynga(1987) cat- 
alogue of open clusters brighter than M v — —4.5.) Our 
value of the CFR corresponds to a starformation rate in 
clusters of 5.9 ± 0.8 10 2 M /Myr within a region of 600 
pc from the Sun for an CIMF with a slope of a — 2 and 
an adopted lower cluster mass limit of WOMq. This cor- 
responds to a surface formation rate of the galactic disk 
near the Sun of 5.2 ± 0.7 10~ 10 Msyr-^c" 2 . 

This value can be compared with the present total star 
formation rate in the galactic disk near the Sun. Lada & 
Lada (2003) derived a SFR of 7 to 10 x 10- 10 M Q y^pc -2 
from embedded clusters in the solar neighbourhood. This 
value is a factor 1.3 to 1.9 higher than the value de- 
rived from the clusters in the Kharchenko sample. The 
difference is most likely due to the fact that many of the 
star clusters formed in embedded clouds will be dispersed 
within 10 Myr. At later ages they would not be recognized 
as clusters. So if the sample of embedded stars studied by 
Lada & Lada (2003) is complete, the infant mortality rate 
of clusters in the solar neighbourhood is about 40 percent. 

We can also estimate the infant mortality rate of the 
clusters from the data in the Kharchenko et al. (2005) 
catalogue. The mean cluster formation rate in the agebin 
of 6.6 < \og(t) < 6.8 is 1.3 10~ 6 clusters/yr, whereas it 
has dropped to 5.4 10~ 7 clusters/yr at 6.8 < log(f) < 7.0. 
This indicates a survival rate of about 40% and an infant 
mortality rate of about 60%. The same rates are found 
if we compare the mean value of the cluster formation 
rates in the age bin of 6.5 < log(f) < 6.9 with that of 
6.9 < log(t) < 7.3. So, the comparison of the cluster for- 
mation rate with the star formation rate of Lada & Lada 
(2003) and the comparison between the formation rates of 
the youngest to the slightly older clusters both suggest an 
infant mortality rate of about 50%. 

6.5. The burst between 250 and 600 Myr ago 

The predictions shown in the top panel of Fig. 1101 and dis- 
cussed above do not explain the bump around log(f/yr) ~ 
8.5 which is higher than any of the distributions for a 
constant cluster formation rate. So the cluster sample of 
Kharchenko et al. (2005) suggests that there was an in- 
creased cluster formation rate around that time. We have 
modelled this with Eq.^lfor several non-constant CFRs. 
The best fit is shown in the lower part of Fig. which 
was calculated for to — 3.3 Myr and with a CFR that is in- 
creased by 0.40 dex between 250 and 600 Myr ago. The fit 
matches the data well. This suggests that the cluster for- 
mation rate was a factor 2.5 higher during this age range. 
Taking into account this burst we find that the mean CFR 
within 600 pc from the Sun during the last Gyr in the so- 
lar neighbourhood was 9.1 ± 3.5 10 2 M Q /Myr which cor- 
responds to a surface formation rate of 8.1 ± 3.0 10 -10 
Mgyr^pc- 2 . 

Zaritzky & Harris (2005) found a peak near 400 Myr 
in the SFR of the SMC. During this peak the SFR was 
at least twice as high as the quiescent SFR of the SMC. 



They attribute this peak to the perigalactic passage of the 
SMC and the Galaxy. Possibly this passage also triggered 
the increased SFR in the galactic disk. 

De la Fuente Marcos & de la Fuente Marcos (2004) 
have studied the starformation history in the solar neigh- 
bourhood, based on various open cluster catalogues. They 
identified five bursts at 0.35, 0.70, 1.13, 1.50 and 1.93 Gyrs 
respectively. The burst that we found between 0.25 and 0.6 
Gyr may correspond to the one of 0.35 or, more likely, to a 
combination of those at 0.35 and 0.70 Gyr found by de la 
Fuente Marcos & de la Fuente Marcos (2004). Our anal- 
ysis of the Kharchenko et al. (2005) cluster sample does 
not confirm the other bursts. 

6.6. The upper mass limit of the clusters 

The theoretical fits of the predicted age distribution to 
the observed one in Fig. ^] shows that the best fit is 
reached if the maximum mass of the clusters formed within 
600 pc from the Sun during the last few Gyr was about 
10 5 Mq. This upper limit nicely explains the steep drop 
in the age distribution at t > 1 Gyr. To check the ro- 
bustness of this conclusion we have also calculated mod- 
els with higher mass upper limits for the clusters. For 
instance, if the upper mass limit was 10 6 Mq then we 
would expect the following predicted and observed num- 
bers of clusters in the oldest logarithmic agebins, indi- 
cated in a vector (logi m i n , logi up ; nr predicted, nr ob- 
served): (8.9,9. 1;5. 1,6), (9.1,9.3; 4.1,3), (9.2,9.4; 3.3,2), 
(9.3,9.5;2.9,1). We see that in these oldest agebins the 
number of observed clusters is within the statistical la un- 
certainty of the number of predicted clusters. This means 
that we cannot exclude the possibility that the maximum 
cluster mass was higher than 10 6 Mq and possibly as high 
as 10 6 Mq. 

An alternative way to consider this point is to find the 
clusters in the Kharchenko sample with the highest initial 
mass. We use the estimate of the present cluster masses, 
derived from the number of member stars, as described in 
Sect. 6.1, and than applied Eq. 0to convert the present 
mass into the initial mass with a disruption parameter 
to = 3.3 Myr. We find that the two initially most massive 
clusters within a distance of 600 pc, i.e. the ones with 
{log(f),log(M)} = {9.2,3.5} and {9.4,2.8} in Fig. had 
an initial mass of 2.5 10 4 and 3.2 10 4 Mq respectively. For 
a CIMF with a = 2 the number of cluster decreases with 
mass as N ~ M~ 2 , so the expected number of clusters 
initially more massive than e.g. 10 5 Mq will be smaller 
than 1. 

The observed mass upper limit is probably determined 
by statistical effects (Hunter et al. 2003; Gieles et al. 
2005b). Its expected value can be estimated from the num- 
ber of observed clusters within 600 pc. For a CIMF with 
a slope of -2 the maximum mass expected in the sample is 
roughly M max ~ N x M; n , where N = 114 is the number 
of observed clusters and M m i n ~ 1O 2 M . So we expect a 
maximum initial mass, set by statistical effects, of about 



H.J.G.L.M. Lamers et al.: Disruption of star clusters in tidal fields 



13 



1 10 4 Mq. (In reality it should be slightly higher because 
the observed number of clusters is already affected by dis- 
ruption.) This agrees quite well with the maximum initial 
mass of about 3 10 4 Mq derived in the previous para- 
graph. 

7. Discussion and summary 

We have derived a simple analytical expression for the 
mass loss from star clusters due to stellar evolution and 
disruption as a function of time, Eq. This expression 
agrees excellently with results of TV-body simulations of 
clusters in the tidal field of our galaxy. The expression 
is derived for mass loss by stellar evolution using the 
GALEV cluster evolution models (Schulz et al. 2002; 
Anders & Fritze-v. Alvensleben 2003) but can easily be 
applied to other cluster evolution models, provided that 
the mass loss due to stellar evolution can be expressed 
by an analytic approximation (e.g. of the type proposed 
in Sect. 3.1). Our analytical expression for the mass loss 
from star clusters is different from the one by Vesperini 

6 Heggie (1997) because they assumed that tidal effects 
decrease the mass of a cluster linearly with time. 

Our method is based on the fact that the disrup- 
tion time of clusters, defined as tdi s = (<iln M/dt) -1 
depends on the mass M, as t d [ s — t (M/ Mq) 1 with 

7 = 0.62 for disruption by two body relaxation in a tidal 
field. This dependence was found both empirically from 
a study of cluster samples in four galaxies by BL03 and 
by Lamers, Gieles & Portegies Zwart (2005), and theo- 
retically from A^-body-simulations by Baumgardt (2001) 
and BM03. The description contains a normalization pa- 
rameter, to, that depends on the environment in the par- 
ent galaxy of the clusters. This parameter can vary by 
more than an order of magnitude between cluster samples 
of different galaxies, as shown by BL03. Simple theoret- 
ical predictions, numerical simulations and empirical de- 
terminations of cluster disruption in a few galaxies showed 
that t ~ C cnv lO~ 47 (p amb /M pc~ 3 )~ ' 5 , where p amh is 
the ambient density in the galaxy at the location of the 
clusters, and C cnv — 300 - 800 Myr (Lamers, Gieles & 
Portegies Zwart, 2005). The value of to may be shorter in 
interacting galaxies, as suggested by the study of the star 
cluster sample of M51 by Gieles et al. (2005a). 

Using our description of mass loss from clusters by stel- 
lar evolution and disruption, we can predict the resulting 
present day mass, M, and age distributions in terms of 
N(M, t)dMdt of the surviving clusters for different cluster 
initial mass functions and different cluster formation his- 
tories. This is given by Eq.^J An integration of N(M,t) 
over mass for different ages gives the expected age distri- 
bution of the surviving clusters. An integration of N(M, i) 
over age for different masses gives the expected mass dis- 
tribution. This method can be applied to any cluster for- 
mation history! 

A comparison between predicted and observed distri- 
butions of selected cluster samples can be used to derive 
the basic properties of that cluster population, such as the 



cluster formation rate, the cluster IMF, and the disrup- 
tion parameter to- To demonstrate this method we have 
applied it to the age distribution of open clusters in the 
solar neighbourhood within 600 pc based on the new clus- 
ter catalogue by Kharchenko et al. (2005). Tests showed 
that the cluster sample in this catalogue is unbiased up 
to a distance of at least 600 pc and possibly 1 kpc. The 
predicted age distribution agrees very well with the em- 
pirical one (Tig. I10[) . The main uncertainty in the derived 
disruption time to is the unknown lower mass detection 
limit of this cluster sample. We estimated this mass limit 
from the given number of cluster members brighter than 
V = 11.5 magn. and found it to be between 80 and 100 
Mq . With these values, the fit of the predicted age distri- 
bution to the observed one shows that the disruption time 
to is 3.3^* Myr, which corresponds to a disruption time 
of a 10 4 Mq cluster of 1.3 ±0.5 Gyr. This is about a factor 

5 shorter than predicted by iV-body simulations of clusters 
in the tidal field of the solar neighbourhood (BM03). The 
difference is possibly due to the fact that BM03 adopted 
a rather high initial central concentration of the clusters 
that is more applicable to globular clusters than to open 
clusters. Moreover encounters with giant molecular clouds 
may also shorten the lifetime of open clusters. 

The present star formation rate in the solar neighbour- 
hood within 600 pc, derived in this paper, is 5.9 ± 0.8 10 2 
M Q /Myr during the last few Gyrs. This corresponds to a 
surface formation rate of 5.2 ± 0.7 10~ 10 M© yr _1 pc~ 2 , 
which is about a factor 0.5 to 0.7 smaller than derived from 
the formation rate of stars in embedded clusters (Lada & 
Lada, 2003). This suggests that a considerable fraction of 
the embedded clusters will be dispersed before reaching 
an age of several Myr. 

The observed age distribution clearly shows evidence 
for a bump in the cluster formation rate between 0.25 and 

6 Gyr, when the formation rate was a factor 2.5 higher 
than before and later. This corresponds to two of the 
bumps in the star formation rate derived from cluster sam- 
ples by de la Fuente Marcos & de la Fuente Marcos (2005). 
The observed bump might be due to an encounter between 
the SMC and the Galaxy (e.g. Zaritsky & Harris, 2005). 

Although the upper mass limit of the observed clusters 
in the solar neighbourhood is less than about 10 4 Mq , we 
show that this does not exclude the possible formation of 
higher mass clusters. Given a cluster IMF with a slope of 
—2, the absence of clusters more massive than 10 4 Mq is 
in agreement with the statistical uncertainty, even if the 
real upper limit for the initial mass was as high as 10 6 

Mq. 
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